Symmetry-breaking Fermi surface deformations 
from central interactions in two dimensions 
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We present a mean field theory of the Pomeranchuk instability in two dimensions, starting from 
a generic central interaction potential described in terms of a few microscopic parameters. For a 
significant range of parameters, the instability is found to be pre-empted by a first-order quantum 
phase transition. We provide the ground state phase diagram in terms of our generic parameters. 



I. INTRODUCTION 

A central theme in the study of strongly cor- 
related electron systems is the appearance of 
novel types of ordering, and phase transitions 
leading to such unconventionally ordered states. 
Phases with liquid crystalline symmetry have 
emerged as an intriguing theme. They have been 
studied and proposed, for example, in the con- 
text of quantum Hall systems, 1-8 s-wave pair- 
ing of polarized fermions, 9 and in Hubbard- 
like models. 10-15 One prominent mechanism for 
such phases is via Pomeranchuk instabilities, 
which are distortion instabilities of the Fermi 
surface. 16 A Pomeranchuk instability occurs in 
the angular-momentum channel I when the corre- 
sponding Fermi liquid parameter Fi is sufficiently 
negative. 16 The phase resulting from an / = 2 
type instability is a nematic, because the orien- 
taion symmetry of the continuum or the orienta- 
tion symmetry of the lattice point group is bro- 
ken, modulo an inversion symmetry, while trans- 
lational symmetry remains unbroken. Pomer- 
anchuk instabilities have received significant at- 
tention recently, both in the continuum 17-23 and 
lattice 15,24-27 contexts. 

In this Article we focus on continuum sys- 
tems, where the Pomeranchuk instability breaks 
a continuous symmetry. This is particularly rel- 
evant to ultracold fermionic gases and low den- 
sity 2D electron systems. The Hamiltonian most 
prominently studied for this case has been of the 
quadrupole-quadrupole type explicitly designed 
in Ref. 17 to produce an I — 2 instability. 21-23 
It is therefore important to study Pomeranchuk 
instabilities arising from more generic Hamilto- 
nians. In a previous paper, two of the present 
authors studied shape deformations of the Fermi 
surface of a three-dimensional system, arising 
from central interactions. 18 In particular, Ref. 18 
finds that a central interaction, if it has a sharp 
feature at a finite length scale ro > kp , can 
cause deformations of the Fermi surface. Non- 



monotonic, "delta-shell" and monotonic, "hard- 
core" repulsive potentials were analysed, and 
they were all found to lead to the effect. A 
screened Coulomb interaction, on the other hand, 
was found not to lead to a Pomeranchuk insta- 
bility. 

In this paper we provide a mean-field treat- 
ment of two dimensions, where much of the cur- 
rent interest lies. In the interest of providing 
generic results, we will parametrize our central 
interaction V{r) by a small number of param- 
eters, namely the values of angular-momentum 
components of the potential and its momentum- 
space derivatives at the Fermi momentum. These 
turn out to be the essential parameters for the 
description of Fermi surface shape distortion 
transtions. A phase diagram in terms of these 
parameters provides much more general informa- 
tion than the consideration of particular forms of 
V(r). 18 Our framework has the added advantage 
that a single description treats not only the tran- 
sition in the I = 2 channel, but in every angular 
momentum channel I > 1. In other words, we in 
fact present a mean field theory not only for dis- 
tortions leading to nematic symmetry, but also 
to the others shown in Fig. 1. 

One important result of our analysis is that 
we find the shape-deformation transitions to be 
of first order for significant regions of parameter 
space. This implies that in many realistic cases 
where one might get shape deformation instabil- 
ities of the Fermi surface, the transition is dis- 
continuous and does not involve quantum critical 
behavior. 

We introduce the model in Sec. II and the mean 
field treatment in Sec. III. The details of the 
theory are worked out in the next three sections. 
Sec. VII describes the parameter regimes where 
we have first-order transitions. In Sec. IX we 
provide a comparison with the Hamiltonian of 
Ref. 17, which is the dominant Hamiltonian used 
in the recent literature for the continuum I = 2 
Pomeranchuk transition. Finally, in Sec. VIII we 



2 




FIG. 1: Shape of Formi surface of 2D continuum 
system, before (thinner line) and after (thicker line) 
a Pomeranchuk instability. Several possible symme- 
tries (I = 2, 3, 4, 5) are shown for the Pomeranchuk 
order. 

use our mean field theory to compute the Landau 
Fermi liquid parameters in terms of the micro- 
scopic interaction potential. In the final section 
we lay out our conclusions. 

II. MODEL 

We consider the following continuum Hamilto- 
nian, describing interacting fermions in two di- 
mensions: 

H = J dripl [e Vr /» - m] 

+ 1 J drj dr'V(\r - r'|)V4v£Vv Vr (1) 

Here e^/ft = p 2 /2m gives the 'bare' disper- 
sion relation in terms of the momentum oper- 
ator p = (ft/i) V r (free fermions of mass m). 
V(r) is a central interaction potential. For sim- 
plicity we consider spinless fermions. (In any 
case, when the interaction is spin-independent 
the Pomeranchuk instability will be degenerate 
in the spin channels. 18 On the other hand for 
spin-dependent interactions there can be emer- 
gent spin-orbit coupling 28 and non-trivial spin 
configurations in momentum space. 29 ) In recip- 
rocal space, Hamiltonian (1) takes the form 

k 

+ ^ 2 T/ ( k " k ') ?/ 'q/2-k^q/2 + k^q/2 + kV' q /2-k, 
k,k',q 

( 

where V fc = / d 2 r fT 1 /^*- 1 ^ and V(k) = 
/ d 2 re~ lkr V(r). The area of the sample is fi. 

There is a large body of literature on behav- 
iors that can emerge from Hamiltonians of the 
type given by Eq. (2). For the attractive case 
see, for example, the discussions of unconven- 
tional superconductivity, Bose-Einstcin conden- 
sation of "preformed pairs" and "fermion con- 
densation" in Refs. 30-33. Here we are con- 
cerned with the emergence of anisotropy of the 



Fermi surface due to repulsive interactions. As 
in the case of anisotropic pairing in superconduc- 
tors, the key is the dependence of the interaction 
strength V(q) on momentum transfer, ftq. In 
the case of a contact potential, V(q) = constant, 
the only possibilities are isotropic: a conventional 
Bardccn-Coopcr-Schrieffer (BCS) instability for 
attractive interactions or a quantum gas-liquid 
transition for repulsion. 

With more general momentum dependence, 
Hamiltonian (2) can lead to an anisotropic state 
through a Pomeranchuk instability. The relevant 
order parameter is 

<&t («/,)) =^ cos (Z0 k ) (4Vk), 

k 

where I — 2,3,4,..., leading to the Fermi sur- 
face shape deformations represented in Fig. 1. 
The I = 1 instability is forbidden in a Galilean- 
invariant system. 18 ' 34 

Pomeranchuk order can be described as 
particle-hole pairing with center-of-mass momen- 
tum q = but finite internal angular momen- 
tum quantum number, I > 0. In this sense, it is 
the analog of "exotic" superconductivity, where 
I > pairing occurs in the particle-particle chan- 
nel. Indeed that analogue has been shown re- 
cently by one of us to extend quantitatively to 
the disorder dependence of a Pomeranchuk tran- 
sition temperature. 35 

It is useful to introduce the partial wave de- 
composition of the interaction potential: 

oo 

V(k - k') = J2 v i( k > k ') cos [I (#k - #k')] (3) 

1=0 

where k = |k|, etc. The amplitudes Vi are 

Vi(k,k') = (1 + M47TX 

/ drr Ji{kr)Ji{k'r)V{r). (4) 
Jo 

(Note that V(k,q) = V{q,k)). We provide be- 
low a mean-field theory that is applicable to any 
Hamiltonian of the form (2) . The theory is quite 
independent of the details of the interaction po- 
tential V(\r — r'|), but depends only on a few pa- 
rameters derived from it. In particular, we will 
use: (i) the values of the amplitudes Vi on the 
Fermi surface (ii) their derivatives on the Fermi 
surface. These are the parameters intro- 
duced in Eq. (21). 
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III. MEAN FIELD THEORY 

Our mean-field theory is based on the following 
ansatz for the ground state: 

|A k ) = n[ ( £ O + e(- £k )^ t (k)]|o). (5) 

k 

Here 0) is the vacuum and V^(k) creates an elec- 
tron in a plane wave state. The ansatz wavefunc- 
tion is a Slater determinant of plane waves. £k is 
an arbitrary dispersion relation. Its form dictates 
which plane wave states are occupied: 

£k = £|k| - M - A k < 0, (6) 

where we have introduced a 'mean field' Ak which 
is the difference between the 'bare' dispersion re- 
lation efe — fi and the renormalized one. 

The mean field Ak is our variational parame- 
ter. This field renormalizes the bare electronic 
dispersion relation, and therefore the Fermi sur- 
face, which is defined as £k = 0. Minimizing the 
ground state energy, 

E= (A k \H |A k ) =^(e k -^)A r k 

k 

+ ^E 7V *[nO)-nk-k')], (7) 

k,k' 

with respect to the functional form of Ak yields 
a self-consistency equation that determines Ak- 
Proceeding as in Ref. 18, we find 

A «= ^£[^(k-q)-^]iVk, (8) 

k 

where V = J d 2 rV(r) and 7V k = (V'kV'k) = 
(— £k)- Note that A q coincides with the 
one-loop, Hartree-Fock approximation to the 
Fcrmionic self-energy S(q, u) (the frequency- 
dependence dropping out at this level, for static 
interactions). When the above equation has more 
than one solution, the one that minimises 




must be chosen. The above expression results 
from substituting Eq. (8) into Eq. (7). The 
second term inside the brackets may be inter- 
preted as a double-counting correction to the 
naive mean-field theory which emerges from the 
variational calculation. 

IV. SMALL FERMI SURFACE 
DEFORMATIONS 

Let us split the mean field Ak into two parts: 
a rotation-symmetric part and a symmetry- 
breaking part. The latter will have a number 



of components, corresponding to different values 
of the angular momentum quantum number of 
the electron-hole pair. Nevertheless near an in- 
stability of the isotropic state, or a sufficiently 
weak first-order transition out of it, we can as- 
sume, save accidental degeneracies, that only one 
of these components is finite. We thus write 

A q = A o (g)+Ai(g)cos(/0 q ), (10) 

or, equivalently, 

£ q = £o(<?) - A ( (<7)cos(^ q ), (II) 

where eo{q) = h 2 q 2 /2m — [i — A (g) is the renor- 
malized dispersion relation before the instability 
sets in and I = 2, 3, 4, . . . determines the symme- 
try of the instability (see Fig. 1). We have chosen 
a particular orientation of the deformation of the 
Fermi surface, without loss of generality. 

In the symmetric phase (zero deformation), the 
Fermi momentum hk F is defined by 

e (k° F ) = 0. 

In the symmetry-broken phase, this quantity de- 
pends on the direction of k, given by the angle 
9: 

k° F -» k F {9) = k° F + 5k F (6). (12) 

The offset of the Fermi vector is given by 

e (k° F + Sk F ) - A z (k F + 8k F ) cos (W) = (13) 

All results presented in the remaining of this pa- 
per have been obtained by solving this equation 
for small deformations of the Fermi surface, i.e. 
under the assumption that 

\Sk F (9)\<^k° F (14) 

in all directions 9. In particular we will assume 
5k F to be small enough that the symmetric part 
of the dispersion relation can be linearised: 

s (k F + 6k F ) KhvpSkp. (15) 

This is quite distinct from the work of other au- 
thors, where non-linear terms in the dispersion 
relation were invoked to stabilize a quantum crit- 
ical point. 17 ' 36 ~ 38 That is discussed in detail in 
Sec. IX. Similarly we will assume that, within 
Sk F of the Fermi vector, the |k|-dependent am- 
plitude of the deformation potential can be ap- 
proximated by a constant: 

A; (k F + 5k F ) « A; (k F ) = A (16) 

Note that we are not writing explicitcly the de- 
pendence of A on the angular momentum quan- 
tum number, I. 
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The approximations (15,16) are valid when the 
following conditions hold for n = 1,2,3,...: 

—L-eWSkZ, ^A {n) Sk^ « hv° F , (17) 

where A< n ) = d n ki{k) / dk n \ k=k o p and e[ n) = 
d n eo(k)/dk n \ k=k o^. Under these conditions, 
Eq. (13) yields the deformation of the Fermi sur- 
face as 



Sk F (9) = A-cos(Z0). 



(18) 



Note that Eq. (18) allows us to re-write the 
small-deformation condition (14) in the form 

A <C hK° F v%. (19) 

Using the decompositions (3) and (10) in the 
self-consistency equation (8), we can express the 
self-consistency equation as 



1 P< 



k° F +Sk F (0) 



dk k Vi(k,k° F ) 
= A. (20) 

In general, this is a self-consistency equation de- 
termining the values of A that minimize (and 
maximize) the energy. The self-concistency 
comes in through the dependence of Sk F (8) on 
A, Eq. (18). Now, the integral equation (20) can 
be reduced to a polynomial equation in A by in- 
voking condition (13) again to keep a finite num- 
ber of terms of the expansion of the interaction 
potential around the Fermi surface: 

Vi(k,k° F ) » V + V'(k-k° F ) 

+ \v"{k-k F f 

+ l i V"'(k-k F ) 3 + ... (21) 

Here V = Vi(k F ,k F ) gives the strength of 
the coupling on the Fermi surface and V' = 
0V5(fc,A$)/0fc| fc=fc o = dV l (k° F ,k)/dk\ k=k o F is its 
slope. Higher derivatives give the curvature, etc. 
Substituting Eq. (21) into Eq. (20) we obtain the 
following equation: 



1 



2tt 



A = — \ Vk° F I d6cos(W)Sk F (6) 



+ 



V + V'k F 



nS) /-27T 



/ d6cos{W)5k F {0f 
Jo 



+ 



y J decos(W)6k F (6) 3 + (22) 

Substituting the value of 8kp(6) given by Eq. (18) 
and carrying out the integration with respect 
to 9 we obtain, to a given order in Sk F /k F = 
A/hk F v F , a polynomial equation in A. 



V. INSTABILITY EQUATION 

The instability equation is found by solving the 
self-consistency equation (20) at lowest order in 
our small parameter expansion, Eq. (21). Using 
Vi(k,k F ) w Vi(k F ,k F ), the critical Vi required 
for the instability is found to be 



V l (k F ,kp) = ^ = V, 



h,p 



crit • 



(23) 



This is our instability equation. Note that 
Vi(kp,k F ) is the only parameter of the interac- 
tion potential entering the instability equation 
(although the I = amplitude is also important, 
as it renormalises the Fermi velocity v F ). For 
Vi(k F ,k F ) < U cr it, the symmetric Fermi surface 
is a (local) energy minimum. For Vi(k F ,k F ) > 
Vent, it is a maximum. 

Writing Vi(k F ,k F ) and v F explicitly in the 
above equation, we find it in the form 

fc2 roc 

-=/ drrV{r)[Mk F r) 2 - Mkprf] (24) 
m Jo 

for I = 1,2,3,.... Note that the above equation 
lacks any solutions with I = 1, as was found in 
D = 3, 18 due to the Galilean invariance of the 
system. 34 

As discussed in Ref. 18 for the three- 
dimensional case, Eq. (23) is a microscopic ver- 
sion of the Pomeranchuk instability condition, 16 
and Eq. (23) reduces to the Pomeranchuk condi- 
tion if we use our mean field theory to compute 
the Landau Fermi liquid parameters in terms 
of the microscopic interaction V(|k — k'|) (see 
Sec. VIII). 



VI. ORDERED STATE 

We now discuss the evolution of the amplitude 
of the deformation, A, in the ordered ground 
state realised when V m (kp,kp) > U C rit- To de- 
scribe the ordered state, we will need to go be- 
yond the lowest order in 5kp/k F in Eq. (21). To 
derive the instability equation (23) , we only used 
the first parameter in the expansion of the Z th 
component of the interaction potential, Eq. (21): 
Vi(k,k F ) V. Let us now also keep the first 
derivative: 



Vi{k,k F )^V + V (k-k° F ). 



(25) 



When substituted in the self-consistency equa- 
tion in the form of Eq. (20), we get, after inte- 
grating with respect to k, Eq. (22). Substituting 
Eq. (18) and carrying out the integrals we get 

U-^A+ — -^3- A 3 . (26) 



Anhvp 



16irh 3 v° F 
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This equation admits two solutions: the trivial 
one, A = 0, minimises the ground state energy 
when V < Krit- On the other hand when V > 
V CT it the above solution is a maximum. If V < 0, 
the minima are at 



A 



hv° F k° F 



± 



iV r , 



V'k F 



1/2 



V - Vci 

Vcr\t 



1/2 



(27) 

The amplitude on the Fermi surface of the defor- 
mation potential thus grows in second-order fash- 
ion, with critical exponent = 1/2, as expected for 
this mean- field theory and depicted in Fig. 2 (a) . 
If V > 0, on the other hand, we obtain the 
unphysical result that the Fermi surface defor- 
mation decreases as the instability point is ap- 
proached from below [see Fig. 2 (b)]. This is 
an indication that linearizing the interaction in 
Eq. (21) is no longer adequate to describe a tran- 
sition which potentially is becoming first-order. 
To address this we go beyond the assumption of 
Eq. (25) in the next section. 



VII. FIRST-ORDER QUANTUM PHASE 
TRANSITION 
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Within our linearized theory we have shown 
that, if V > 0, unphysical solutions occur. These 
solutions are suggestive of the isotropic state be- 
ing a local minimum with the true ground state 
separated from it by a first order transition where 
the Fermi surface jumps to one of lower symme- 
try. Linearisation of Vi(k,k F ) around k = k F , 
Eq. (25), is therefore not realistic for change that 
is no longer infinitesimal near the transition. We 
therefore carry out the expansion (21) to the 
third, rather than first, order. It is useful to in- 
troduce the following, dimensionless parameters 
describing the interaction potential and the re- 
sulting deformation of the band structure: 



(n) . V K F 



V'crit 

A 

hv° F k° F 



(28) 
(29) 



In terms of these, the self-consistency equation 
takes the form 



1 / v"\ v'" 
(v~l)\+-(vi+ V Y )\3 + V —\^0. (30) 



Taking v" 



and solving for A recovers 



Eq. (27) as A = ±2 (-l/i/) 1/2 (v - 1) 1/2 . How- 
ever, Eq. (30) has valid solutions for v' > 0, too. 
All the solutions are straightforward to obtain 
analytically. Two examples are plotted in Fig. 2. 



FIG. 2: Values of the amplitude on the Fermi sur- 
face of the deformation potential, A, that minimise 
the free energy (solid lines), plotted as functions of 
the coupling strength in units of its "critical" value, 
VVVcrit. (a) Vkp/Vctit = -1 and V {n) k n F /V cvit « 
for all n = 2,3,... [Eq. (27)] (b) V'k F /V clit = 1, 
V 2 k 2 F /V C rit = -2 + 1/50, V 3 fcf7Krit = -1 and 
V {n) k^/V C rit w for n > 4 [Eq. (30)]. The dashed 
line in panel (b) indicates an additional stationary 
point, but it is a maximum, not a minimum. The re- 
gion where the transition takes place has been blown 
out in the inset. The dotted lines indicate the criti- 
cal coupling of Eq. (23) [panels (a) and (b)] and the 
lower bound in Eq. (32) [panel (b)]. The structure of 
the free energy in the different parameter regions has 
been sketched for illustration. 



We note that there are two types of behaviour. 
For 



v' + v"/2 < 



(31) 



(which reduces to v' < in the limit of small 
v") there is a second-order transition as we found 
above. However, when the above condition is not 
met, there is a range of values of V for which 
the free energy has a triple-well structure. The 
instability is then pre-empted by a weakly first- 
order quantum phase transition: a small "jump" 
in the shape of the Fermi surface. The value 
of the coupling strength at which this happens, 
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FIG. 3: Phase diagram for small, symmetry-breaking 
deformations of the Fermi surface with a given value 
of the angular momentum quantum number I. In 
the symmetric state, the Fermi surface is circular. In 
the broken-symmetry state, it has one of the config- 
urations of Fig. 1. In the shaded region one of the 
two states is stable, and the other meta-stable, so 
the transition takes place through a first order jump. 
The parameter V" k%/V CI it = —1/20, which controls 
the width of this region [Eq. (32)]. 



^jump 



< V CI it, is bounded by 



3(v' + v"/2Y 
4v'" 



< 



^jump 
Vr.rit 



< 1. 



(32) 



The corresponding phase diagram is shown in 
Fig. 3. The exact location of the jump de- 
pends on additional parameters characterising 
the fermion-fermion interaction (Specifically, it 
depends on the off-diagonal second and third 
derivatives of Vi(k,q) with respect to k and q.) 

The above is only valid for v'" < 0. If v'" > 0, 
higher-order terms describing the dependence of 
Vi(k,q) on k and q even further away from the 
Fermi surface become important. In that case 
the transition is no longer even weakly first-order 
and the precise behaviour of the model depends 
on more specific details. Naturally as the whole 
solution relies on the assumption that 8kp{9) is 
small the range in which it is reliable is restricted 
to values of the parameters for which the jump is 
small. 

It is important to stress that the above ar- 
guments rely on the approximation of Eq. (18). 
In the present case, however, the approximation 
made in that equation requires not only that 
| A' | < v F , but indeed that Eq. (17) holds for 
n = 1,2,3. Equivalcntly we have neglected the 
additional renormalisation of the Fermi velocity 
and higher-order derivatives of the dispersion re- 
lation resulting from the symmetry-breaking part 
of the Fermi surface deformation, compared to 
the symmetric contribution. 



VIII. LANDAU PARAMETERS 

Variation of Eq. (7) gives the change of the 
ground state energy associated with an arbi- 
trary re-distribution of the fermions in momen- 
tum space, iVk — > iVk + SN^: 



SE = Y^ ekSN k 



1 

2TL 



^/(k,k')<W k cW k ,. 

k,k' 

(33) 



This coincides with the usual expression from 
Landau Fermi liquid theory. 44 The effective 
interaction between quasiparticles is given by 
/(k,k') = V - V(k-k') in terms of the mi- 
croscopic parameters. The Landau Fermi liquid 
parameters can be defined from this function in 
the usual way: 

oo 

f(k,k') = — ) Y,Flcos[l(6 k -9v)} 

for k=k' = k° F (34) 

where p(0) = l/2irhv F is the density of states at 
the Fermi energy. Thus 



Ft = p(0) (S L0 V - Vi) 



(35) 



In terms of these Landau parameters, the Pomer- 
anchuk instability equation (23) of section V 
takes the usual form, Fi < —2. The microscopic 
parameters V, V" , etc. introduced in Sec. VI are 
related to derivatives of /(k, k') at the Fermi en- 
ergy and can thus be regarded as generalisations 
of the Landau parameters. 



IX. QUADRUPOLE-QUADRUPOLE 
INTERACTIONS 



In Ref. 17 Oganesyan, Kivelson and Frad- 
kin (OKF) introduced an effective Hamiltonian 
which has been employed widely to study the 
/ = 2 instability on the continuum. 22,23 It fea- 
tures an anisotropic, "quadrupole-quadrupole" 
effective interaction. By contrast our Hamilto- 
nian, Eq. (1), features a central interaction po- 
tential. In this Section we relate the two. 

The OKF Hamiltonian is 



Hokf 
1 



h 2 v 2 



ip(r) 



2m 

dr I dr'^ 2 (r-r')tr[Qt(r)Q(r')]. (36) 



This corresponds to postulating directly an 
anisotropic density-density interaction with I = 2 
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symmetry which, by analogy with classical liq- 
uid crystals, is assumed to be of quadrupole- 
quadrupole form. 39 The "quadrupoles" used here 
are the quadrupole moments of the electronic mo- 
mentum distribution, 17 




</>(!•) . 



The philosophy behind the effective Hamilto- 
nian in Eq. (36) is that in order to describe 
the important fluctuations on approaching the 
Pomcranchuk instability, in a given channel I, 
from the isotropic state it is not important to 
include fluctuations tending to create either dis- 
tortions of the Fermi surface shape with different 



symmetry or an inhomogeneous state (backscat- 
tering terms). To understand the relationship of 
theories based on this approximation to our anal- 
ysis of the generic central Hamiltonian of Eq. (2), 
we will construct a central Hamiltonian whose 
form is constrained so that its leading instability 
has an effective Hamiltonian of the OKF form. 

In order to compare with the generic Hamilto- 
nian (2), we first move to reciprocal space: 

Q(r) = £ V t (ki)Q ka ^(k 2 )e i < kl - ka >- r . 



kik 2 



One finds Ql 1 = 



-Ql 2 



-(k 2 /k F ) cos(26), 
and Q 12 = Q 21 = (k 2 /k 2 F ) sin(20), so that 
tr[Q kl Qk 2 ] = 2 (k 2 k 2 /k F ) 0082(0! - 9 2 ). Thus 



k b k,k',q 

COS [2 (6> k _ q/2 - 6 , k'+q/ 2 )] ^ k+q/2 ^_ q/2 ^ k -q/2^k'+q/2 (37) 

where the Fourier transform of the interaction potential is defined by ^(r — r') = 
fi_1 E q eiq (r_r ' ) ^2(q) and we have introduced the notation 

?k = £k + 7^ ]T T2 (<0 ( k - k " cos I 2 (^- q - 0k)] • (38) 
F q 

Like our Hamiltonian, Eq. (37) features a pairwise interaction that preserves the total momentum 
of the pair, q. Note, however, the complicated dependence on q, k' and k. By contrast in Eq. (2) 
the interaction depends only on the transferred momentum, h(k — k'). The dependence on q, in 
particular, means that the OKF interaction is not uniform. However the actual form of .7-2 (q) is not 
very important 17 and it is customary 17,26 to take ^ 2 (q) = ^2^,0- With this assumption Eq. (37) 
takes the simpler form 

Hokf = ^>~k - /MV-k - ^J-2k 2 k' 2 cos[2(0 k -0 k O]V k 4'^k' (39) 

k F k,k' 

I 



Let us now compare this Hamiltonian to the form 
considered here, given by Eq. (2). In our trial 
ground state, Eq. (5), the only terms in the triple 
sum that contribute to the energy have k = k' or 
k = — k'. Omitting all other terms from Eq. (2) 
it takes the simpler form 

H = Y ( £ k - M) V'kV'k 

k 

+ i E t y ( k ' - k ) - ^l^' • (4°) 

k,k' 

Substituting in this expression the partial wave 
expansion of the interaction potential, Eq. (3), 



we obtain a scries of different interaction terms, 
labelled by I, which lead to Fermi surface de- 
formations with different symmetries. Near a 
second- or weakly first-order Pomeranchuk dis- 
tortion with I = 2, only the corresponding term 
need be considered (see above). Neglecting all 
the others the Hamiltonian takes the form 

H = Y ( £ k - M) V'lV'k 

k 

+^ E V2 ^ k ') cos [ 2 ( 0k - 0k ')i vM'V>kW, 
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whose interaction part coincides with that in 
Eq. (39) if we take 

V t (k,k') = -^T^' 2 - (42) 

Thus although the OKF Hamiltonian docs not 
correspond to a central interaction potential, 
within our theory it would give the same re- 
sults as a hypothetical central interaction, whose 
I = 2 component happens to be given by Eq. (42). 
Note that inserting Eq. (35) into Eq. (42) we ob- 
tain T-2. oc i<2, as expected. 17 

From this result we note that, since Ti is neg- 
ative at the instability, it follows from Eq. (42) 
that 

V = -M^>0, (43) 

Kp 

(and that V" > too). This implies that within 
our mean field theory, based on a linearised dis- 
persion relation, Eq. (15), the central potential 
model which captures the OKF Hamiltonian is in 
the parameter regime where the Pomeranchuk in- 
stability is actually first order. Indeed to stabilise 
a quantum critical point for the OKF Hamilto- 
nian it is essential to include non-linear terms 
of the symmetric dispersion relation, £o(k), as 
noted in Rcfs. 17,36-38. Conversely, our comple- 
mentary approach shows the generic conditions 
under which a quantum critical point can be sta- 
bilised without invoking such non-linearities of 
the dispersion relation. This is achieved instead 
by properly taking into account ultra-violet cut- 
offs implicit in any given central interaction po- 
tential, V(r). 

Our results are consistent with other work con- 
sidering specific microscopic realizations of Fermi 
surface instabilities which also found wide regions 
where the transition is first order 15,26 . We stress 
that the main lesson one should extract from this 
is that the order of the quantum phase transi- 
tion is a very delicate issue, depending on fine 
details of the effective interaction and the band 
structure. Indeed the higher-order terms of the 
dispersion relation alluded to above modify the 
coefficients of A 3 and A 5 in Eq. (30) [for a de- 
tailed analysis, see Ref. 36.] For example, a large 
enough 6k p term in the expansion of So(k F +5kF ) 
can change the sign of the coefficient of A 3 , and 
hence the order of the transition. Thus a full 
analysis going beyond present calculations would 
have to treat non-linearities in the interaction 
and the dispersion relation on an equal footing. 40 
Moreover higher-order effects beyond the scope of 
mean field theories may well upset this balance 
one way or the other. Such non-mean field ef- 
fects will certainly become important whenever 



the condition (19) does not hold — e.g. if the 
phase transition is of first order, but not weakly 
so. 



X. CONCLUSION 

We have provided a mean-field theory for con- 
tinuum Pomeranchuk transitions in two dimen- 
sions. The theory is expressed in terms of a few 
pertinent parameters (V, V , V" , ...) for each 
angular momentum channel. This makes the the- 
ory quite general, and applicable for a wide class 
of central interactions in which the symmetry 
breaking is not put in explicitly by hand. 

Our main results are Eqs. (31) and (32), 
which determine the phase diagram in Fig. 3. 
They apply to any central interaction poten- 
tial in a two-dimensional continuum, described 
in terms of a few dimensionless parameters, de- 
fined in Eq. (28). Depending on the form of 
the interaction, our theory may lead either to a 
first- or second-order quantum phase transition. 
Thus our approach is complementary to other 
work 17 ' 36-38 where a quantum critical point was 
stabilised by non-linear terms in the dispersion 
relation. 

A continuum theory is useful for several rea- 
sons. One is the direct relevance to several exper- 
imental systems where a Pomeranchuk transition 
might be realized, for example: 2D electron lay- 
ers at semiconductor hetero junctions or on liq- 
uid helium; layered helium systems 41 (where a 
new phase intervening between the Fermi liquid 
and Mott insulator states has been observed 42 ); 
or in a cold-atom setting with trapped fermionic 
atoms with dipolar repulsion. These are situa- 
tions where underlying lattice structures are not 
expected to play a role. Theoretically, continuum 
Pomeranchuk transitions are fascinating because 
the resulting broken symmetries are remarkable. 
While some effects of the nematic symmetry (bro- 
ken 0(2)/ Z2) have already been explored, 17 we 
believe there are further implications, for exam- 
ple, effects of nematic half-vortex excitations in 
such a medium. 43 Our framework also provides 
a description of transitions in higher-/ channels 
(leading to 0(2) /Zi symmetry broken states), 
which presumably leads to a broader class of 
interesting excitations and properties. Some of 
these issues are currently under investigation. 
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